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Abstract: Based on minimal supersymmetric standard model, neutralino 
dark matter annihilation induced gamma ray emission from galaxy cluster 
1E0657-56 is calculated. The merge of bullet-like subcluster with the main 
cluster is also investigated. 
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1 Introduction 

Both astronomical and cosmological observations have concluded that Dark Matter(DM) 
contributes a quarter of the energy partition in universe and DM dominates the matter 
composition in galaxiy. Among various theoretical DM particles candidates, neutralino 
predicted by SUSY models is the most attractive one. Though there is no electromag- 
netic radiations, as the antipaticle of itself, they can annihilate into a pair of photons 
or other particles which eventually decay to the final state photons. Theoretically, the 
corresponding gamma ray flux is proportional to the square of the dark matter density. 
Therefore, such a signal can most probably be observed from the center of the galaxies 

or substructure in galaxies. 

Galaxy cluster 1E0657-56 is a very massive cluster who has an exceptional bullet-like 

subcluster in it. The weak lensing observation of bullet cluster indicates the separation 
of baryonic mass component and unseen matter induced by the merger, which provided 
a direct empirical proof of the existence of DM[lJ. This encourage us to study the cor- 
responding DM signals from it. In this work, we calculate the DM annihilation induced 
7-ray flux from the central direction of main cluster and bullet cluster, respectively. 



2 Gamma Ray Flux from Dark Matter Annihilation 

Dark matter particle is usuallay considered as the Light Supersymmetric Particle(LSP) 
which, in most supersymmetric scenarios, is regarded as the neutralino x- Theoret- 
ically, gamma ray from \X annihilation is expected to exhibit both monochromatic 
and continuum spectrum while the later part makes the dominant contribution and 
therefore becomes the major concern of this work. Suppose that the dark matter is 
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concentrated in a spherical halo of virial radius r v { r and denote the density profile p x 
as a function of the radial distance, then energy spectrum of gamma ray from DM 
annihilation in a galaxy can be written as[2j 
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where Nj(E) = Y2f describes the differential energy spectrum of the gamma ray 

from DM annihilation for all final states F, d refers to the distance to the galactic 

center. Af2 represents the solid angle for a given angular resolution of the detector. 

Actually, Eq.([I]) can be divided by two parts, one called "particle factor" which is 

exclusively determined by the feature of particle, and the other called "cosmological 

factor" which decided by the density distribution of the local dark matter halo. This 

allow us to calculate them separately. 

In the Minimal Super symmetric Standard Model(MSSM), certain SUSY parameter 

space gives the information of the mass and the interaction strength of DM particles. 
Generally, SUSY parameter space is decided by seven independent SUSY parameters 
p, M2, Ma, rrig, tan/3, A t and A/,. For the particle factor part, we perform a random 
scanning of these parameters within the corresponding ranges by using the package 
DarkSUSY[3j. The value of these parameters are constrained by the theoretical consis- 
tency. Furthermore, the result is checked by current accelerator data and experimental 
observations. One of strong limits comes from the observational data of the relic abun- 
dance of the cold dark matter, VtcDhih 2 = 0.113^Q QQg, specifically [3J. We make use of 
3cr error bar range to limit the parameter space, i.e., 0.086 < Q x h 2 < 0.137. The relic 
bundance of neutralino particle less than the minimal limit indicates a subdominant 
DM component. Then, a rescaling of DM density will be done as p(r) — > £p(r) with 

ti = nh 2 /(n x h 2 ) mm ^. 

3 Gamma Ray Emission from Galaxy Cluster 1E0657-56 

The main cluster of galaxy cluster 1E0657-56 is found that the halo density contribution 
is well fitted by NFW[6] profile with concentration C200 = 3.0 and virial radius r2oo = 
1.64/i _1 Mpc, corresponding to a virial mass and velocity of M200 = 1-5 x 10 15 M Q and 
^200 = 2380A;ms _1 [7]. The bullet substructure, located about 0.48/i _1 Mpc to the 
centre of the main cluster, is the most massive substructure in the cluster with the 
M200 = 1.5 x 1O 14 M , C200 = 7.2 and r 200 = 0.78/i _1 Mpc. The NFW density profile 
can be described as 

Px{r) = (r/r s )(l+r/r s )2' (2) 
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where p s and r s are the scale density and scale radius respectively. These two free 
parameters can be determined by the measurements of the virial mass and the concen- 
tration parameter of the halo. In order to avoid the unphysical singularity when r — > 0, 

we take a core radius r m j n = 10 -11 kpc, within which the DM density keep invariant [8j . 
Cluster 1E0657-56 is a very massive cluster formed by a plenty of sub-clusters. 

Dark matter halo clustered in small length scales was wildly discussed |H1 EH [HI IT2] . 

Since clumpy sub-halo will greatly enhance the 7-ray emission, the effect including 

substructure should be considered. The clumpy magnitude of dark halo is mainly 

determined by the product of two parameters, / and 8. where / is defined as the fraction 

of clumpy component in the total dark halo, and 5 is a dimensionless parameter which 

gives the effective contrast between the sub-halo density and the local halo density po, 

8 = Pcl/po ~ 0(1000). The factor / is estimated to be around between 5% ~ 20% in 

the literature |13 |, 115]. Following the approach of Ref. [TT] and the references therein, 

the cosmological factor {J(iP))aci can be rewrite as 

where we take the mass fraction of sub-structure as f ma in = 10% and f bullet = 50% [16] 
in our calculation. 

4 Results and Discussions 

Based on the theoretical frame given above, dark matter induced 7-ray flux from centre 
of main cluster and bullet subcluster are calculated, respectively. The central distance 
between the two cluster is about 0.48/i _1 Mpc which is equivalent to the angle of 
6.322 x 10~ 4 between the light of sight. Furthermore, the overlap of the two cluster 
is considered where the angular resolution of AS7 = 10~ 3 is applied. For the main 
cluster, we get (J(VO)ac = (3.61 x lO" 3 )^^ + (6.56 x 10~ 4 ) overlap- While we get 
(J"(V0)Afi = (1-56 x 10~ 3 )dumpy + (2-79 x W~ 3 ) ove rla P for the bullet subcluster. This 
shows that the two have almost the same intesity of gamma ray emission. We consider 
the constrains come from both particle physics and astrophysics and the final result is 

presented in Figure [TJ 

The result shows that, for the neutralino annihilation induced 7-ray emission, 

galaxy cluster 1E0657-56 has a double peak structure on the spacial distribution. 
For those Imaging Air Cherenkov Telescopes(IACTs), their sensitivity can achieve 
~ W~ 12 cm^ 2 s^ 1 in the lOOGeV ~ WTeV region. Even though, there are still sev- 
eral orders of magnitude difference to be detectable, which is hard to reach simply by 

increasing statistical time. 

In comparison with cherenkov technique, EAS experiments such as MilagrofTT] and 

ARGQ[l8] have the advantage in large field of view and higher duty cycle which enable 

them to posses the capability in observing simultaneously a large number of sources of 

sub-structures. Therefore, they have the potential in improving the detection sensitivity 

of gamma ray from DM annihilation. 
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Figure 1: The integrated gamma ray flux from the centre of galaxy cluster 1E0657-56 
above the threshold energy of 1,10, 100, 500, lOOOGeV, where NFW halo profile and the 
angular resolution of AQ = 10~ 3 is set. Each neutralino mass refers to a set of SUSY 
parameter space which checked by all the theoretical and experimental limit. 
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